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In a previous paper, we presented solution to the Slavnov-Taylor identity for the QCD effective 
action, and argued that the action terms containing (anti)ghost fields are unique. These terms have 
the same form as those in the classical action, but the gluon and (anti)ghost effective fields are 
convoluted with gluon and ghost dressing functions Ga and G c , the latter containing perturbative 
and nonperturbative effects (but not including the soliton-like vacuum effects). In the present 
work we show how the perturbative QCD (pQCD) can be incorporated into the framework of 
this action, and we present explicit one- loop pQCD expressions for Ga and G c . We then go on to 
check the consistency of the obtained results by considering an antighost Dyson-Schwinger equation 
(DSE) . By solving the relations that result from the Legendre transformation leading to the effective 
action, we obtain the effective fields as power expansions of sources. We check explicitly that the 
aforementioned one-loop functions Ga and G c fulfil the antighost DSE at the linear source level. We 
further explicitly check that these one-loop Ga and G c have the regularization-scale and momentum 
dependence consistent with the antighost DSE at the quadratic source level. These checks suggest 
that the the effective action with dressing functions represents a consistent framework for treating 
QCD, at least at the one- loop level. 

PACS numbers: 11.15.Tk, 12.38.Aw,12.38.-t,12.38.Bx 



I. INTRODUCTION 



Nonperturbative methods are essential in order to rigorously prove well accepted QCD phenomena such as con- 
finement and chiral symmetry breaking. In principle the effective action contains all these effects, allowing for the 
calculation of all proper (one-particle-irreducible) n-point Green functions. Furthermore, its minimum gives the true 
quantum mechanical vacuum of the theory. One way to obtain information about the effective action is through the 
Slavnov-Taylor (ST) identity |||, §, which is a consequence of the invariance of the classical (tree level) action 
with respect to the BRST transformation ^, |j. The ST identity for the effective action relates different functional 
derivatives of this action with respect to the effective fields. In the N = 1 four-dimensional supersymmetric theory 
without matter, the form of the solution to the ST has been found in Ref. Q. In Ref. 0, we presented a solution to 
the ST identity for the effective action of nonsupersymmetric non-Abelian gauge theory. 

QCD is an example of nonsupersymmetric non-Abelian gauge theory. The solution has the form of the classical 
(QCD) action, but the fields entering it are convolutions of the corresponding effective fields with pertaining dressing 
functions. These functions are supposed to contain the quantum effects of the theory - perturbative and nonpertur- 
bative, with the possible exception of soliton-like vacuum effects. Further, in that work we argued that the presented 
terms of the solution which contain the ghost (and antighost) fields are the only ones possible. Stated differently, 
the terms in the effective action containing (anti)ghost terms have the classical form and all the quantum effects 
connected with Green functions with at least one (anti)ghost leg must be contained in the ghost and gluon dressing 
functions Ga and G c . 

We also assume that those terms in the effective action solution which contain no (anti)ghost fields have essentially 
the classical form, but with the fields convoluted with the dressing functions. Further, the coupling parameter is 
running. In this case, we have not presented arguments that such a solution to the ST identity for the effective 
action is the only possible one. The main attractive feature of this effective action is that all the quantum effects are 
parametrized with a finite number of dressing functions and with the running function of the coupling. 

A natural question appearing is whether the known QCD results can be consistently incorporated into the afore- 
mentioned effective action, i.e., into the gluon, ghost and quark dressing functions. In this work we show how the 
known one-loop perturbative QCD (pQCD) results are reflected in this effective action. In fact, we find explicit 
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one-loop pQCD expressions for the ghost and gluon dressing functions Ga and G c by requiring that the the correct 
known one-loop gluon and ghost propagators be reproduced starting from our effective action. It is not clear in 
advance whether our QCD effective action framework with the obtained dressing functions gives a self-consistent 
description of QCD in the high-momentum (perturbative) regime. We address this question of self-consistency in the 
following way. We first write down an antighost Dyson-Schwinger equation (DSE) for QCD. Then we express the 
effective action fields as a power series in the (external) sources, by solving the relations that follow from the Legendre 
transformation leading from the path integral to the effective action. These power series involve the aforementioned 
dressing functions. We then insert these expansions in the antighost DSE, with the one-loop dressing functions Ga 
and G c . We explicitly check that the regularization-scale and the momentum dependence of these Ga and G c are 
consistent with the antighost DSE at the linear source and the quadratic source level, 1 while at the linear source level 
we even check the consistency of the constant ("finite") parts of the one-loop Ga and G c . All this strongly suggests 
that our effective action framework is consistent with pQCD in the corresponding (high) momentum region. It further 
keeps open the realistic possibility that the dressing functions of our effective action represent a sufficient framework 
for including not just all the (high momentum) perturbative, but also (lower momentum) nonperturbative quantum 
effects. 

In Sec. [lj we recapitulate the results of Ref. M and present the aforementioned antighost DSE. In Sec. Ill 



we 



deduce the one-loop gluon and ghost dressing functions Ga and G c from requiring the correct one-loop gluon and 



ghost propagators. In Sec. [V we write down the relations arising from the Legendre transformation leading to the 
effective action, and obtain the effective fields as power expansions in sources. In Sec. [v] we explicitly check that the 
obtained one-loop dressing functions Ga and G c are consistent with the antighost DSE at the linear level of sources. 
At the quadratic level of sources, we explicitly check that the regularization-scale and momentum dependence of the 



one-loop Ga and G c is consistent with the antighost DSE. In Sec. VI we summarize our results 



II. THE EFFECTIVE ACTION AND THE ANTIGHOST DYSON-SCHWINGER EQUATION 

In order to fix notations, we first write the well-known classical action of QCD, in the Lorentz gauge 
S QC b[A, b, c, q, q] = j d 4 x { - ^ Tr [F^(A(x))F^ (A(x))} + q {j) i ^V^A)q U) 

~ Tr {[d^(x)] 2 ) - 2 Tr [ i b(x)d^ V„(A) c(x)} } . (1) 

Here, A, b, c, and qj are the gluon, antighost, ghost, and quark (j'th flavor) fields, and a is the gauge parameter. 
For simplicity we omit the quark masses. For all the fields X in the adjoint representation (X = A, b, c) we use the 
notation X = X a T a , where T a — X a /2 are the gauge group generators (N c x N c hermitian matrices in the color 
space) with the normalization 

Tr (T a T b ) = i(5 afc , (T a ) t = T a , (2a) 
[T b ,T c ] = if abc T a . (2b) 

The antighost and ghost fields are hermitian Grassmann numbers w = b, = c. The covariant derivatives in the 
ghost part in (Q) are 

V„(A)c = d^c + i[A^,c] , (3a) 

= d„A v - d v A lL + i[A M , A,} , (3b) 

^^(A)q(j) = dpqy) + iA^q^) . (3c) 

The fields X a , q^, the coupling g, and the gauge-fixing parameter a in the action (|l|) are bare, i.e., the theory has a 
large (but finite) UV momentum cutoff A 

X^X^, q (j) =q$, 5^ 5(A), a = a™ = 5 2 (A)^ , (4) 
and the cutoff is implemented through a gauge invariant regularization, e.g. dimensional regularization. 



The two DSE's involve indirectly three- and four-point connected Green functions. 
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A solution to the Slavnov-Taylor (ST) identity for the effective action T has a similar form as the classical action 
(|I|). However, the fields are modified 1^1 |. 2 In particular, the terms containing the (anti)ghost fields and the 
quark fields have just the same form as in the classical action, but with JhJ, and the pure gauge boson terms 
have gauge invariant form, but with A i— » A 



T QC n[A,b,c,q,q] = J d 4 x [cf (A(x)) + Cf (q U) (x), A{x)) 

-i Tr {[d^{x)f) - 2 Tr [ i b(x) 9^V M (i) c(x)] } . (5) 

Here, C\ s (A(x)) is a gauge invariant combination of the gluon field A. The modified fields X appearing in (||) are the 
original fields convoluted with the corresponding dressing functions Gx 

■ - i ( n ... / ,7 J ,.' ,-• i- , 



A;(x) = (G-/ o A°) (*) = J d*x' G A \x - x')^(x') , (6a) 

c. a {x) = (G- 1 oc a )(x)= J d 4 x' G^ 1 (x — x')c a (x') , (6b) 

9C0«( X ) = (Gq 1 °1u)) s (x) ^ J d A x' G^ix-x'^q^ix') , (6c) 

6 a (x) = (G A o6 a )(a;)= / dVG A (i-i')6°(i') . (6d) 



Notice that the dressing functions Gx appearing in relations (6a) and ( |6d| ) are the same (Ga)- In ([k]), the spinor 
structure is reflected in indices = 1,2,3,4. The spacetime structure of the dressing functions is (x\Gx \x') — 
Gx(x — x') due to translational invariance. These functions are real for X = A, c, and are in general complex 4x4 
spinor matrices for X = quy Furthermore, the Fourier transforms of the dressing functions Ga and G c 

G x (k 2 ) = J d A x exp(ik-x) G x {x) (7) 

are functions of k 2 due to Lorentz invariance. As a consequence, Gx{x) — Gx(-i) (X — A,c). Apart from the 
aforementioned invariance properites of the dressing functions, the ST- identity for the effective action does not tell 
us anything about these functions. 

We have argued in Ref. |7j that the terms containing the (anti)ghost fields in the solution ([3]) of the ST identity have a 
unique form as given in Eq. (||) . This was a consequence of the assumption of Rcf. 0] that the term ~ Lc 2 of the effective 
action fulfils separately the ST identity for the effective action. Here, L is an auxiliary (nonpropagating) background 
field which couples at the tree level to the BRST transformation of the ghost field c in the form 2Tr[L(x)c(x) 2 ] in the 
(classical) Lagrangian density. The aforementioned assumption was then shown to imply the following structure of 
the ~ Lc 2 term of the effective action: 

d 4 x 2 Tr [(G c o L) (x) (G" 1 o c) (x) (G" 1 o c) (x)] , (8) 

where G c (x) is, a priori, an arbitrary dressing function. An argument was given in Ref. Q in favor of the afore- 
mentioned assumption. Furthermore, arguments were presented there which suggest the structure (||) as the only 
consistent ~ Lc 2 structure even when the mentioned assumption is not adopted. The structure of this correlator was 
then shown in Ref. Q to restrict the gluon-ghost-antighost term to the one written in Eq. (^), and the physical part 
of the effective action was shown to be gauge-invariant in terms of the effective fields convoluted with the dressing 
functions (|) X = G x l o X. 

The gauge-fixing term in (^) is also unique, due to its known insensitivity to quantum corrections. In (^), the 
expression £| which contains the quark fields is gauge-invariant combination of q^, q^ and A. In our further 
considerations in this paper, the precise form of these terms will not matter as we will concentrate on the question of 
consistency of the (anti)ghost and gluon sectors of the effective action (|^) in the perturbative region. 

The pure gluon term C\ (A(x)) in the effective action has a form somewhat different from the form of the pure 
gluon term in the classical action (fil). This should be expected because the latter term is the only one which contains, 



2 The exception is the gauge-fixing part which is unchanged under quantum corrections. 



4 



in addition to the field(s), the coupling parameter g. This parameter is running. Thus, it is natural to expect that 
the latter parameter in the effective action is not just normalized to a specific value g(po), where po is a specific, 
in principle arbitrary, renormalization scale. Rather, when going into the Fourier-transformed (momentum) space, 
we should associate to each fc-mode of the dressed gluon field A the value of g(k). Specifically, while the pure gluon 
sector in the classical action (Q) has the following form in the momentum space: 



S#[A] 



1 



d 4 k 



Tr [F^ u (A)(k) F^(A)(-k)] 



(9) 



where F } u(A)(k) denotes the Fourier transform of F^(A)(x) 
action (0) contains the analogous term 



F^(A(x)), the pure gluon sector in the effective 



[A] = -- 



d 4 k 



(27T)4 ff 2(_ fc2) 



Tr 



F, u {A){k) F^(A)(-k) 



(10) 



However, in contrast to (Q), the expression (|T^) is not gauge invariant when A t— > A, although the full gluon sector 
of the effective action (HJ must have this property as a consequence of the ST-identity. Thus, gluon contributions 
additional to those in ( |l0| ) must appear. In order to illustrate how to obtain these additional gluon terms, we will 
find them explicitly in the one-loop perturbative approximation. In this case, we have the relation 



1 



1 



A) 



9 2 {-k 2 ) 
where 0o 



1 (lliVc 



In 



-k' 2 



12 



4tt 2 
- 2n/) . 



Mo 



(11a) 
(lib) 



Here, /3q is the one-loop coefficient for the beta- function for a s = g 2 / (4 ir), N c = 3 is the number of colors, n/ is the 
number of active quark flavors. Using the one-loop evolution form (11a) in (|Io|), we can rewrite the latter in the form 
of the following spacetime integrals: 



r 



1— loop 

gU 



[A] = 



2 4tt 2 



d A x Tr 
d^x Tr 



F^A)(x) F^(A)(x) 



<fi 



F MV (A)(x) In ) F^(A)(x) 



(12) 



The integrand in the second term is not gauge invariant when A i— > A. However, the minimal extension of the latter 
integrand to recover the aforementioned property is to replace <9 M i— > V M (-A) 



t1— loop 



[^4] 



1 



2ff 2 (A*o) 

1 / ,il 

2 4tt 2 



d*x Tr 



F MV (A){x) F^(A){x) 



d x Tr 



F^(A)(x) In 



v 2 (A(x)) 



Mo 



F» v (A)(x) 



(13) 



By expanding the new logarithm in ( |l3| ) around its value at A — 0, we can immediately see that the new terms, i.e., 
the difference of the integrands of the second integrals in (|l2|) [«-> (fu^)] and (|l3|), are terms at least cubic in gluon 
fields and of the type 1/9, i.e., in momentum space these are power-suppressed terms of the type 1/k. Such terms 
will not be relevant in our perturbative considerations, only the terms that depend logarithmically on the momenta 
will be. Therefore, in the perturbative considerations, we can consider the form (|l^) to represent the entire gluonic 
sector of the effective action. 

In Appendix |A| we show explicitly how Eq. ([l2]) follows from Eqs. (|l0|) and (|ll]), present clear expressions for 
the logarithmic operators appearing in the integrals in Eqs. ( |l2| ) and (|l3|), and show (tilde-)gauge invariance of the 
expression ([l3]). Further, in that Appendix we argue how to proceed in the general (nonperturbative) case, i.e., when 
l/g(—k 2 ) is a general function of /c 2 , in order to obtain from ( |To| ) the (tilde-)gauge invariant version of the gluonic 
effective action. 

For further use in the manuscript, we will write here a general Dyson-Schwinger-type equation (DSE) which the 
QCD path integral has to fulfill. The QCD path integral Z = exp(— iW) is defined as 



exp(-iW[J,p,ri,j,j)) = [dA][db][dc][dq(j)][dq U )] exp (iS Q cT>[A,b, c, q, q] + iE[A, b, c, q, q; J, p, 77, j, j]) , (14a) 
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with 



SL4,6, c,q,q;J,p,r),j,j] 



2 Tr 



d 4 x J^A^ 1 + i / d A x p{x)b{x) + i / cTa; r)(x)c(x) 



d 4 x (q U) j U) +3{j)q{j)) , 



(14b) 



where J, p, 77, j, j are the (external) sources for the corresponding fields. The path integral Z ([14a]) is just a function 
of sources, with the fields integrated out. Therefore, in particular, Z is invariant under the (infinitesimal) shift 
of the antighost fields b(x) <— * b(x) + e(x). This invariance can be rewritten, after some straightforward algebraic 
manipulations, 3 in the form 



d 2 c a {x) - f abc d" 



A»(x)c c (x) + 



S 2 W 



SJ^(x) 5r] c {3 



-p a {x) . 



(15) 



Moreover, due to the relation 8W/8rj c {x) = —ic°(x) which follows from the Legendre transformation that connects 
W and Tqcrj, this relation can be rewritten as 



d 2 c a (x) - f abc 



-P a (x) 



(16) 



Here it is understood that the effective fields c and A are functions of sources. We will call the identity (|T^) the 
antighost Dyson-Schwinger equation (&-DSE). The usual 6-DSE equation is obtained by applying to (16) the variation 
5/(Sp d (y)) and then setting all the sources and fields equal to zero. In (fl6|), the fields and the sources are in general 
nonzero. 



III. PQCD PREDICTIONS FOR THE DRESSING FUNCTIONS 

We will here deduce some implications of perturbative QCD (pQCD) for the hitherto unknown dressing functions 
appearing in (^J). We will write explicit expressions for them as predicted by one-loop pQCD. 
The pQCD prediction for the full gluon propagator is 



5 2 W 



(A:) 



5 2 T 



(k)=g 2 (A) Z 3 (A,K) D«l<$ K Xk), 



where Z 3 (A,K) appears as a factor in the gluon field rescaling function for the theory with cutoff K (K 2 = —k 2 

A^=Z g (K,K)Zl /2 {kJ<) Af) , 
and Z g (A, K) is the rescaling ("running") function for the QCD coupling parameter 

g(A) = Z g (A,K) g{K) . 



(17) 



(18) 



(19) 



Here we use a simplified notation for the running g: g(K) = g(—k 2 ), g(A) = g(A 2 ). The rescaling functions Z3 and Z g 
are calculable in pQCD. In fll7|), (fc) is the usual gluon propagator in the theory with cutoff K (K 2 = — k 2 ), 



i.e., the usual tree level propagator 



fcjll kfi 2 j:(K) fyfl ^jJ2 

k 2 5 k 2 



(20) 



and ^ K > = a( K '/g 2 (K) (~ g°) is the usual gau ge p arameter. On the other hand, the gluon propagator as predicted 
by the effective action (||), with the gluon part (|10|) , is 



S 2 T 



&Ap y SAfj^ 



(k) = 5 aic 



(-1) 
k 2 



h h 



g 2 (K) (G A (k 2 )) ( Wa - ) + 



h h 

i /. 1 rv 11 



k 2 



(21) 



Care should be taken about the Grassmannian anticommuting character of p and b. 



6 



We stress here that this is the full gluon propagator in our framework. The additional nonperturbative terms in the 
gluonic sector of the effective action which appeared in Eq. jl3| ) as a consequence of the substitution d i— ► V(A(x)) in 
the logarithm, are terms at least cubic in the gluon fields (<~ (A(k)) 3 /k) and thus do not affect the gluon propagator 
which is determined solely by the ~ (A(k)) 2 terms of the gluon part (|To|). Comparing the transversal parts of the 
pQCD propagator jl7|) and of our propagator (El]), we conclude 



G A (k 2 ) = (j^j Z 1 3 /2 (A,K) = Z g (A,K)Z 1 3 /2 (A,K). 



(22) 



The longitudinal parts of the pQCD propagator (|T^) and of our propagator (pl| ) agree now automatically, because of 
the relation (^) for the gauge parameter a and the known relation 

^=Z 3 (A,K) ^ . 

Now we analogously consider the ghost propagator. The pQCD prediction for the full ghost propagator is 

5 2 W \ ( 5 2 T 



\Srj ai 6p a2 J v ' \5b a2 Sc a i J y ' 12 2 v ' k 2 



(23) 
(24) 



where (A, K) is the ghost and antighost field rescaling function for the theory with cutoff K (K 2 = —k 2 

c (A) = Z (c)l/2 (A) K) C (K) t b (A) = Z f X l\A, K) ftW . 

On the other hand, our effective action (||) predicts 

s 2 r . (G c {k 2 y 



\hj — $a 1 a 2 



\G A (k 2 ) 



i 

k 2 



b a2 Sc ai 

Comparing this with the pQCD result (|24|) gives us 

G c (k 2 ) = zi c) (A,K) G A (k 2 ) 

= Z { 2 C) (A,K) Z 3 1/2 (A,^) Z g {A,K) 

Although we will not use the explicit form of the quark dressing function (matrix) G qi we mention here that completely 
analogous considerations involving the quark propagator give us a possible (but not unique) pQCD form for G q 

G q {k 2 ) 6e =zl /2 (A,K) exp(i8(k 2 ))6 6 



(25) 
(26) 

(27) 



Oft/ 



where 9(k 2 ) is a momentum-dependent phase, and Z 2 (A,K) is the quark field rescaling function (K 2 



Y (A) = z i/2 



'Li) 



(28) 
(29) 



The rescaling ("running") functions Z 3 , Z 2 , Z2 and Z g are all computable in pQCD, and at the one- loop level are 
given by 



Z 3 (A,K) = 

Z { 2 C) (A,K) = 

Z 2 (A,K) = 

Z g (A,K) = 



1 + -a s 

7T 



1 + ~a s 

7T 



1 + -a s 

7T 



1 + -a s 

7T 









+ K 3 



+ 0(a;) , 



.„ 1 £1 



.(c)- 



n 2 m — j 



K 2 J 

K 2 , 



K'2 



0(a 
+ 0(a 2 s ) , 
+ 0(a. 



2 ) 

s I 1 



2 ) 

s I ' 



where a s = g 2 /(Ai:) and the coefficients Kx are well known 



.(c) _ 



K'2 



i [(13 - 30N C - 4n f ] 



1 (TV 2 - 1) 



2iV c 



T/3o 



24 



(HJV c -2n/) 



(30a) 
(30b) 
(30c) 
(30d) 

(31a) 
(31b) 
(31c) 
(31d) 
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Here, ./V c = 3 is the number of colors, n/ is the number of active quark flavors, and (3$ is from Eq. (lit). Further, the 
"finite" one-loop parts, which are A- and fc-independent, are known as well. For example, in the MS renormahzation 
scheme, using dimensional regularization, we have 



«3 



f-(i-f) + jd-«) 2 



Nr. 



10 



(32a) 
(32b) 



The constants £2 and k g are also known, but we will not need them at the level of our analysis. The "finite" 
part constant (32a) is written for the case when the active quarks are massless. The two constants ( |32| ) can be 
obtained by evaluating the dimcnsionally regularized gluon and ghost self-energy parts - see, for example, the book 
H [Eqs. (A. 21) and (A. 23) there]. In the above formulas, we denoted by A the MS UV cutoff which is obtained by 
the formal substitution [2/(4— D) — 7 + ln(47r)] 1— > ln(A 2 ), where D (— > 4) is the dimension in the regularization. The 
coupling a s is, in principle, the "bare" MS coupling a s = a s (A 2 ;MS) here. However, the momentum dependence of 
a s in expressions ( pC| ) is irrelevant at the ~a s level as it affects only ~a 2 terms. 

Inserting the known one- loop results (|30|)-(^) into ( p^ ) and (p7f), we obtain the one- loop pQCD predictions for the 
gluon and ghost dressing functions G a and G c 



G x (k 2 ) 



1 

Z7T 



Kx In 



A 2 



-k 2 



+ KX 



+ 0(a 2 s ) 



(X = A,c) 



K A = («3-i3o) = --(3 + e)iV c , 

k c = (24 c) + k 3 - A)) = ~^N C , 
KA = K 3 + 2k„ , 



~:N r . 



(33a) 
(33b) 

(33c) 
(33d) 

(33e) 



The freedom of choice of the UV regularization scale A is reflected in the dressing functions Gx which depend on 



A 



G x {k 2 ) = G^\k 2 ) 



(34) 



Up until now, we considered the renormahzation scale in the dressing functions to be a very high UV cutoff scale 
A ^> Aqcd- But the arbitrariness of A implies that, in the Wilsonian sense, we could have started from the classical 
action with a lower cutoff /i < A, where fi is a renormahzation scale. Then, in all our formulas, we would have to 
make the simple substitution A 1— * ji. The resulting rescaling of the dressing functions under the change A 1 — s- is 
then 



A 



Z g (A,») Z 3 1/2 (A, M ) Gf , 
Z^(K^) Zl /2 (K^) Z g (A,ri GM 



(35a) 
(35b) 



It is then straightforward to see that the field A and the combination 6c, appearin g in the effective action (|5|), are 
independent of the choice of the renormahzation scale fi. This can be seen via Eqs. (p5|), ( |l8|) and (|2^), in the latter 
two replacing K 1— > fi. Analogously, it can be shown that the combination Qij)Q(j) appearing in the effective action is 
independent of \i. 



IV. FIELDS AS POWER EXPANSIONS IN SOURCES 



In this Section we will write dow n the transformations associated with the Legendre transformation which relates 
the path integral function W (14a) with the effective action Tqcd @)- We will then solve them, thus obtaining the 
effective fields c and A as power expansions in sources. We recall that the following Legendre transformation must 
relate W with Tqcd: 



T[A, b, c,q,q] 



■W[J,p,ri,j,j] - E[A,b,c,q,q; J,p,rj,j,j] 



(36) 



8 



where S is given in Eq. (14b). From here follows the 6(antighost)-relation 

ST 



Sb a (x) 



St}b ST t \ \ ■ at ^ 
- — o -=-{x) = ip (x) 



By Eqs. (|d|) and (||) we have 



Sb b (y) 
Sb a (x) 



$abG A (y - x) , 



ST 
Sb b (y) 



(>V„(I) ~c(y) 



(37) 



(38) 



Applying to d37j ) the convolution i G A 4 o from the left, and using the aforementioned symmetry G y4 1 (— x) = G A 1 (x), 
we then obtain 



d z c a (x) - f^d" [A" (x)c c (x) 



d 4 x'G~ A 1 (x-x')p a (x') . 



(39) 



This is the Legendre-related b- relation in spacetime coordinates. Applying the Fourier transformation of type (^) to 
this relation, we obtain the following Legendre-related &(antighost)-relation in four-momentum coordinates, involving 
the corresponding Fourier transformed functions: 



-k 2 ^ 



G c (fc 2 ) 



i f abc k» 



d 4 h A h (ki) c c (fc-fci) 



P a (k) 



(2tt) 4 G A (k 2 ) G c ((k - kx) 2 ) ~ G A {k 2 )' 
Completely analogously, we obtain the following Legendre-related c(ghost)-relation in four-momentum space: 

eabc f d 4 kx G A {(k-hf) 



k 1 G A (k ) b a (k)+i f a 



-(k - kxTA b (kx)b c (k - k x ) = G c (k 2 ) v a (k) 



(2tt) 4 G A (k 2 ) 

Very analogous, but algebraically more involved manipulations of the Legendre-related A(gluon)-relation 

ST 



SA°{(x) 



SAf; 2 ST \ 



SAl\ SA a2 



(40) 



(41) 



(42) 
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lead to the following explicit form of the Legendre-related A(gluon)-relation in four-momentum coordinates: 



1 



g 2 (K)G A (k 2 ) \ 
+i f aia2a3 G A (k 2 ) 



d 4 h 



k 2 



9 ^ K) tr 



K\(k) 



[2it) 4 GA(fc 2 )GA((fc-fci) 2 ) 
x [(fcf A a3 ^{kx) - k% 2 A a3fJ,1 (ki)) A%{k - ki) - k^AlKk^A^ik - fci)] 



+f aia2b f a3a * b G A {k 2 ) 



d 4 hd 4 k 2 



(2tt)8 G A (k 2 )G A (k 2 )G A ({k -k x - k 2 ) 2 ) 



A a2 ^{kx)A a3 ^{k 2 )A a Mk -kx- k 2 ) 



d 4 kx 



d 4 kx G A (kl) 



-k^b a2 (kx)c a3 {k - kx) 



(2tt) 4 G c {{k-kx) 2 ) 
(g- - ) (fc 1 ) 7 G- i 1 (fc)7°) T" 1 ^ 1 (G-\k - kx)q(k - kx)) = - G a (k 2 )J^(k) . 



(43) 



The g-relation from the Legendre transformation would involve, analogously, terms ~ q/fi and ~ Aq^ on the left-hand 
side and the quark current on the right-hand side. In the above relations, and in the rest of the article, we omit 
for simplicity any notational reference to the dependence of the effective fields, the dressing functions, and the gauge 
parameter a on the UV cutoff A. 

The Legendre-related relations (40), (41), and (f43| ) connect the effective fields with the sources. They allow us to 
obtain the effective fields as expansions in powers of sources 



- y(i) 



X = X 



(44) 
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Here we denote by the part of the field X of power n in the sources (J£, p, fl,3(j)iju})- We assume that the fields 
go to zero when the sources go to zero {X^ — 0), i.e., we assume that there are no soliton-like vacuum effects. When 
inserting these expansions in the Legendre-related relations, and requiring that each coefficient in the expansion fulfill 
the equation, we obtain an infinite series of Legendre-related equations for X^ n ' } s which can be solved successively. 
For example, the ^-relation (M) thus gives the series of equations 



k 2 



G c (k 2 ) 

k 2 
G c (k 2 ) 

k 2 
G c (k 2 ) 



G A (k 2 ; 



c Q (fc) (2) +i f abc k fi 
c a (k)^ +i f abc k^ 



P a (k) , 
d 4 fci 



1 



^u(fci) (1) c c (fc-^i) (1) = , 



G A (kl)G c ((k - kxY) 



(2tt) 4 G A {kl)G c {{k-k x fY^ 
d 4 fci 

A b {kx) {x) c c {k-kx) {2) + A b u (kx)^c c (k - 



etc. Solving these equations successively for c' 1 ', c^ 2 \ c' 3 ', . . ., we obtain 



c a (fc)« 


G A (k 2 )k 2P (k) ' 




c a (k)W 


= l f abc G c (k 2 )^ J 


r d 4 kx 

(2tt) 4 


c a (fc) (3) 




r d 4 k x 

(2tt) 4 



G A {k\)G c {{k - h) 2 )(k - h)- 



r A b Jkx)^p c (k-kx) , 



G A (k()G c ((k - kx) 2 ) 



A b Jkx) {1) c c {k-kx) {2) + A b Jkx)^c c (k - fer)« 



(45a) 
(45b) 

(45c) 

(46a) 
(46b) 

(46c) 



etc. In an analogous, though algebraically more involved way, by inserting expansions ( [44| ) into the Legendre-related 
^4-relation (0) we obtain expressions for A^> 



2 1 



A^(k)^ = g 2 (K){G A (k 2 )Y P^(k)J^{k) , 



k 2 ' 

GA(k 2 )„ m f d A ki G A (k 2 )G A ((k — ki) 2 ) 



(47a) 



(2tt) 4 k 2 (k-kx) 2 



x [{k^P^Ah) - (fci + fe) M3 P^ 4 (fc 1 )) P% 5 (k - kx)] J^(h)J^(k - kx) + 0( VP ,j um) ) , (47b) 
etc., where we denoted by P^ LV the function proportional to the gluon propagator ( ^ij ) 



9m 



k k 



1 fc^ x fc^ 2 



5 2 (if) (Ga(F)) 2 fc 2 



A- A- 

'I'll 



3mu 



Ml "7*2 ) _|_ ^(if) "7*1 "7*2 



A- 2 



(48) 



In the analysis of the next Section, we will not need the explicit terms @(vP->j(j)j(j)) m (^b)- They originate from 

terms ~ c^M 1 ) and ffyjfffjv We can obtain the coefficient functions b^ n \ cf n \ in a similar way, but we will not 

need them in the following analy si s. Th e cub ic- in -s ourc es part of c can be obtained by inserting in the expression 
(pel) the explicit solutions (|46a|)-(|46bD and (pa|)-(pb|). 



V. CHECKING THE ANTIGHOST DS EQUATION 

In this Section we will perform a consistency check of the effective action framework (||) for the high momenta 
(pQCD) region. In Section [II] wc deduced the dressing functions for the region of high momenta by requiring that 
they, in the framework of (p), should give us the known pQCD two-point Green functions. In particular, the known 
one-loop pQCD expressions for the gluon and ghost propagators resulted in the high-momentum behavior ( |33| ) for the 
dressing functions G A and G c . On the other hand, in the expansion in sources various Dyson-Schwinger equations 
(DSE's) can give us (infinite) series of equations for the dressing functions. For example, the antighost DSE (6-DSE) 
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( |l6| ) gives, after inserting into it for the effective fields the power expansions in sources discussed in the previous 
Section, one equation for each power of sources 



6J£{x) 



= 



(49a) 
(49b) 



etc. In p art A of this Section, we will check explicitly that the aforementioned one-loop pQCD solutions ( p3| ) fulfil 
Eq. (49a). In part B we will check explicitly that their dependence on the UV cutoff A and on the momentum k is 
compatible with Eq. fl49b| ). 

We first apply the Fourier transformations to these relations, obtaining relations for the corresponding Fourier- 
transformed functions 



-k 2 c a (k) w + f abc k" 



(k) = -p a (k) 



-k 2 c a (k)W+ir b w j i^L4(fc 1 ) (1) c c (fc-fc 1 )( i )+r fcc fc^^J(fc) = o 



(50a) 
(50b) 



etc. The variational derivatives of the type (Sc^ /SJ)(k) in Eqs. ( f50| ) denote simply the Fourier transforms, with 
respect to x, of the functions {5c(x)^ /8J(x)). 

We note that the term Sc^/SJ in the linear- in-sources DSE (49a) contains implicitly the three-point connected 
Green fun ction S 2 / (5 pS J) — i S 3 W/(Sp5ri5J) ~ (bcA). Analogously, the term Sc^/8J in the quadratic-in-sources 
DSE (49b) contains implicitly the four-point connected Green function 5 4 W/(Sp57]5JSJ) ~ (bcAA). 



A. Linear-in-sources DS equation 



First we will show, in our effective action fram ewor k (pf) and (Eof) , that th e on e-loop dressing f uncti ons 



consi stent with the linear-in-sources DS equation (50a). When we use in Eq. ( J50a| ) the expressions (46a), ( |46b| ), and 
(47a), the common factor p a {k) factorizes out and we end up with the equivalent relation (we multiply with —G A (k 2 )) 



1 



( A ) d 4 fcx G c {k 2 )G A {{k-k 1 ) 2 ) i 



G c (k 2 ) +ig 2 (K)N c ^ J ^^r h) r ^ P^(k-h)k^=G A (k 2 ) 



(51) 



All the squared momenta in the denominators, here and in the rest of the article, are understood to have -Me added to 
ensure causality. The integral must be understood as having been regularized with a gauge-invariant cutoff A, which 
can be implemented with dimensional regularization (see Sec. HI and Appendix |B|). 

Further, we insert in ( |5l| ) the one-loop results (|33|), and thus are allowed to keep only terms up to order 0(a\). 
This means that we replace in the integral G c {k 2 )GA{{k — &i) 2 ) •— > 1 and ^( K - K i) |_> £( A ) an d obtain the equivalent 
(one-loop) relation 



G A (k 2 )-G c {k 2 ) = i4ira s (K)N c 



k 2 



(A) 



d 4 k! 1 

(2tt) 4 kjik-h) 



rc(A) i^ jk — k^^jk-ki)^ 



fcf .(52) 



I n Ap pendix we explicitly evaluated the integral on the right-hand side of DSE (|52j), with the result given in Eq. 
(B15). On the other hand, the one-loop expression results ( p3| ) for the dressing functions give for the left-hand side 
of DSE (152 



G A (k 2 )~G c (k 2 ) = -a s (K 2 )N c 



1 

167T 



(3-£ (A) )m 



A 2 



0(a 2 s 



(53) 



This is just the same as the result (B15) for the right-hand side of DSE (|52|), at the one-loop level at which we 
are working. This shows explicitly, in our effective action framework, that the one-loop dressing functions j33| ) are 
consiste nt w it h the DSE (|l]); or equivalently with the first (linear in sources) DSE (50a) in the series of the antighost 
DSE's (pa|), (fall), etc. 
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B. Quadratic-in-sources DS equation 



We now proceed to show, in our effective action framework (||) and (|h]), that the A 2 - and k 2 -dependent parts 
of the one- loop dressing functions Gx ( |33| ) are consistent with the quadratic-in-sources antighost DSE (50b). The 
explicit form of the first two terms in this DSE, in term s of the d ressing functions and sources, can be obtained by 
straightforward insertion of the expressions (46a), (46b) and (47a) 



k 2 c a (fc) (2) + i f abc k" 



(A) 



d 4 fci 
(2tt)*"'' 



-i 9 2 {K) k, 



(A) 



d 4 fci 



(2tt) 4 k 2 (k - hi) 



G c (k 2 ) 
G c ((fc-fc 1 ) 2 )G A (fc 2 ) 



- 1 



(/ afcc ^(fciV c (^-fci)) 



(54) 



We recall that in our convention, A ~ g, J ~ g , p ~ g . But the expression in brackets in the integral is ~ g , thus 
the expression (|54|) is ~ g 3 ~ a^ 2 . Using the one-loop expressions ( j33| ) for the dressing functions, the expression in 
brackets is 



G c (k 2 ) 



G c {{k - fci) 2 )G A (fc?) 



- 1 



hi 



A 2 



finite 



0(a 2 



(55) 



Thus, using the explicit form (33b ) for a a, expression 



can be rewritten as 



k 2 c a (fc) (2) +i f abc k 



be 



^(3 + ^)N c a 2 s (-k 2 ) k. 



(2tt) 4 " 



4(fc x )Wc c (fc 
1 



(2tt) 4 fc 2 (fc-fci) 2 



In 



A 2 



-fc 2 



finite 



(f abc Jt(ki)p c (k-h))+0(a 5 s / 2 ) (56) 



The calculation of the last term in DSE ( |50b| ), i.e., its reduction to a form of the type (|56|), is more involved since 
it entails evaluation of (a variational derivative of) the cubic-in-sources fu ncti on c' 3 '. The calculation is given in 
Appendix in the Feynman gauge (^ A ^ = 1), with the result given in Eq. (|C8|). Since the sum in the first brackets 

for the first two terms in 



there is one, we see that the third term result ( |Cq ) is just the negative of the result 
the Feynman gauge. 4 This explicitly shows that the A 2 - and fc 2 -dependence of the one- 



consistent with the quadratic-in-sources antighost DSE (50b) in our effective action framework 



oop dressing functions (33) is 



VI. SUMMARY 

Based on arguments presented in our previous work Ref. 0, we propose a specific form (||) of the QCD effective 
action which basically repeats the structure of the classical QCD action, but where all the fields X are replaced by 
their convolutions with dressing functions Gx'- X = G x x o X. The proposed pure gluon part is given in Eq. jic|), 
and we argue that it should be supplemented by terms which restore gauge invariance in the dressed gluon field A 
and that these supplementary pure gluon terms do not contribute to perturbative effects. In general, the proposed 
effective action should have all perturbative and nonperturbative effects (except for possible soliton-like vacuum 
effects) encoded in the dressing functions Gx- 

We then investigated the consistency of the proposed framework in the regime of application of perturbative 
QCD (pQCD). First we deduced the gluon and ghost dressing functions Ga and G c from the requirement that the 
full gluon and antighost-ghost propagators of the framework should agree with those known from pQCD (with an 
UV regularization scale A), and obtained explicit pQCD expressions for Ga and G c at the one-loop level. Stated 
differently, we fixed the two dressing functions of the framework by requiring the agreement of two two-point Green 
functions calculated in this framework with those of pQCD. But the two dressing functions appear also in many 
other (three-point, four-point, ...) Green functions. These Green functions can be related with the two-point ones via 
Dyson-Schwinger equations (DSE's). A priori, it is far from clear whether a very limited number of dressing functions 
fulfills many (infinite number of) DSE's. We performed explicit checks of the antighost DSE (6-DSE) at the linear- 



4 The "finite", i.e. A-independent parts, in addition to ln(A 2 /(— k 2 )), may also appear in Eqs. (C8) and (j56|), the algebra of including 
them becomes very involved, and they are not considered in the present analysis. 
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and at the quadratic-in-sources level, with the aforementioned one- loop dressing functions G a and G c . These two 
DSE's involve implicitly the connected three-point (bcA) and four-point (bcAA) Green functions. Thus we presented 
nontrivial checks of the consistency of the proposed QCD framework in the perturbative (high momentum) region. 

However, in the quadratic-in-sources DSE we considered only the large logarithms containing the regularization 
scale [ln(A 2 /(— k 2 )] in the dressing functions and in the DSE's, leaving the algebraically involved problem of the 
cancelation of the A-independent ("finite") terms in the DSE's for future investigation. Subsequently, we plan to 
investigate, in the proposed QCD effective action framework, the behavior of the quark sector in the high momentum 
region, and subsequently the behavior of all sectors in the low momentum region. 
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APPENDIX A: LOGARITHM OF THE (COVARIANT) D'ALAMBERTIAN 

In this Appendix we write down explicit expressions for the logarithm of d'alambertian d 2 = 9 M <9 M and the logarithm 
of its (tilde- ) covariant version V 2 (A) = V M (A)V A '(A), and clarify why expression ( |l2] ) is not gauge invariant and 
expression ("J]) is. 

For any function B(x) in spacetime, the function [ln(<9 2 //i 2 )i?](2) can be written in terms of the Taylor expansion 
of the logarithm around the point ln(^io/A*o) = 1 



In 



B 



(x) 



1 



1 



lid 2 -ti)~^(d 2 -ti) 2 + ^{d 2 -ti? + 
o z Mo °Mo 



B(x) 



(Al) 



The Fourier transformation (^) of this function is obtained by integration by parts which replaces d 2 



-k 2 



Mo 



(ft) 



+ Mo) " A(fc 2 + MoT ~ A(fc 2 + MoT 
z Mo Mo 



Mo 



B(k) 



In 



Mo 



B(k) 



(A2) 



This shows that the one-loop relation ( [Tl"|) inserted in expression ( |io|) gives expression jl2]). Further, since d 2 F tJ -"(A)(x) 
is not gauge c ovar iant [F^ v (A)(x) is gauge covariant], expression Tr[F fJil/ (A)(x)d 2 F fJ -"(A)(x)] is not gauge invariant, 
and then Eq. (Al) implies that the second integrand on the right-hand side of Eq. (|lj) is not gauge invariant. 

Unlike d 2 F^ u (A)(x), expression V 2 {A)F^ '(A)(x) is gauge covariant. Therefore, Tt[FuJA)(x)V 2 (A)F^(A)(x)} is 
gauge invariant. The logarithmic version can be written in complete analogy with Eq. ( |Al| ) 



In 



V 2 (i) 



Mo 



B 



-I(V 2 (i) - nl) -^(V 2 (i) 2 - fi 2 ) 2 + ^(V 2 (i) - m 2 ) 3 + 
Mo z Mo J Mo 



B(x) . (A3) 



From here we see that, when B(x) = F tll/ (A)(x), this expression is also gauge covariant, and thus the trace representing 
the second integrand on the right-hand side of Eq. (|l3| ) is gauge invariant. 

From this procedure, we can also see how to proceed in the general (nonperturbative) case, i.e., when l/g(—k 2 ) = 
Q(—k 2 ) is a general function of k 2 , in order to obtain from ( [To| ) the (tilde-)gauge invariant version 



^± c ? (_ fc 2 ) Xr [^(i)( fc ) F^(A)(-k) 



Q(-k 2 



g(-k 2 ) 



d 4 x Tr 



F, 



UA)(x)g(v 2 (A(x))) F^(A)(a 



(A4a) 
(A4b) 



where G{—k 2 ) should be represented as a power expansion. 
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APPENDIX B: REGULARIZED INTEGRALS 



In this article, we need the regularized form of the following types of integrals: 



1 



(2tt)4 kffa+qtf 



/•(A) rf4k 

{l;/ii;Aii/i 2 ;MiM2M3} _ / 1 



(2tt) 4 ^(fc 1+(?1 ) 2 (fci+92) 2 



f1 . 1,^1 . jLMl £.A«2 . 1,^1 ,L,M2 jLM3 1 



(Bla) 
(Bib) 



where the various options are indicated in the curly brackets. The factors in the denominators are understood to have 
+i £ added to ensure causality. The UV regularization is denoted in (Bl) by (A) above the integration symbol. It 
is implemented by dimensional regularization, i.e., in D dimensions, thus replacing e? 4 fci/(27r) 4 i— > dP k\ / {2n) D and 
9^9^ =4h*D. The UV cutoff A is introduced in the final results, after D — > 4, according to the substitution 



4- D 



7 + ln(47r) = ln(A 2 ) 



(B2) 



representing thus the effective "bare" UV cutoff in MS scheme. The relevant formulas leading to the evaluation of 
integrals (Bl) can be found, for example, in Refs. || [|. 

The product of denominators is transformed into one denominator via the method of Feynman parameters 



1 



AxA 2 ...A r , 



Specifically, 



fc 2 (fci + gi) 2 



dx\dx 2 5{x\ + x 2 — 1)- 



£ = k\ + x 2 q\ 



A^=x 2 {l-x 2 ){-q\) ; 



and 



fc 2 (fc 1 +g 1 )2(fc 1 +( Z2 )2 



l r i r i 



JO JO 



dx\dx 2 dx 3 5(xi + x 2 + x 3 — 1)- 



(P- AO))" 

I = ki + x 2 qi + x 3 q 2 , = -XiX 2 qf ~ xix 3 q^ ~ x 2 x 3 (qi - q 2 ) 2 

Introducing the integration variable £ according to ( B4b| ), (B5b), the integrals (Bl) can be rewritten as 

d D £ 1 



(B3) 

(B4a) 
(B4b) 

(B5a) 
(B5b) 



lim 

D^4 



Jo 



x{l;£^-x 2 q^} 

r l pi pi 



dx\dx 2 S(xi + x 2 — 1) 
ah 



(2ir) D (p_ A (2))< 



^(9i,92) {1;Mi;AllAl2;MlM2M3} = 2 Urn / f [ d Xl dx 2 dx 3 5(xi + x 2 + x 3 - 1) 

-^{l;q^;q^q^;q^q^q^}\ , 

> J L ' ^ ^ ' * * ^ ' Iq=l — X 2 qi—X3q2 



d D £ 



(2n) D (P-A(3)y 



(B6a) 



(B6b) 



The integration over £ can be performed after Wick-rotation, using dimensional regularization, and introducing 
e ee (4-I>)/2 



d D £ 




1 


% 


(2ir) D 


(£ 2 


-A) 2 


(4tt) 2 


dP£ 




e 


i 


(2n) D 


(£ 2 


- A) 3 


(4tt) 2 


d D £ 




1 


i 



+ ln(47r) - ln(A) ) + 0(e) 



(2tt) d (£2 _ A )S 



i^(l- 1) ri- 7 + ln(47r)-ln(A)J+0(e) 



(4tt) 2 2 A 



--+0(e) (finite) 



(B7a) 
(B7b) 
(B7c) 
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Further, by symmetry considerations the following identities hold: 

d D £ {l^-J^H^H^} 



(2tt) d (P - A) r 
dPt l^i^ 
(2Vp (I 2 - A)" = D 



, 



d 4 £ 



(2tt) 4 [P - A) r 



(B8a) 
(B8b) 



Applying the identities (|B7|)-(B8) to (|B6|), and at the end performing the limit e = (A—D)/2 — » 0, and then introducing 
the cutoff notation (B2), we obtain the regularized integrals (Bl) 



m 

i{ q y 



(4tt) 2 



In 



2(4tt) 2 



(4tt) 
i 



-q 2 
In 

dx 3 



+ 2 



A 2 



l-x 3 



dxo 



i f 1 



dx 3 



1—2:3 



(4tt) 2 4 



A 2 



F{qi,q2]X2,Xz) 
(x 2 qi+x 3 q%) 
F(qi,q2;x 2 ,x 3 ) 

- -5 M1M2 / dx- 



(B9a) 
(B9b) 
(B9c) 
(B9d) 



da; 2 In 



F(qi,q 2 ; x 2 , x 3 ) 



dx 3 



1 3:3 ^ (x 2 qi+x 3 q 2 ) fll (x2qi+x 3 q2)> 12 
F(qi,q 2 ;x2,x 3 ) 



(B9e) 



/(5i,<? 2 )™ 3 = 



+ - 



(4tt) 2 
1-X3 



- ^ [ff MlMa (9i+?2r +ff M2M3 (91 +9 2 r +fl" 3Ml (9i+92r] In 



dx 2 [g^^qi+x^f 3 + g^^qi+x^Y 1 + g^^qi+x^Y 2 ] In 
(x 2 gi +x 3 g 2 ) Ml {x 2 qi +x 3 q 2 Y 2 (x 2 qi+x 3 q 2 Y 3 \ 



F{qi,q 2 ;x 2 ,x 3 ) 



dxo 



F(qi,q2;x 2 ,x 3 ) 



where we denoted 



F(qi,q 2 ;x2,x 3 ) = -X 2 (l-X2)ql - X 3 {l-X 3 )q\ + 2x 2 x 3 (qi ■ q 2 ) 



(B9f) 



(BIO) 



and t h e squ ared scale (— q 2 ) appearing in Eqs. ( |B9e| )-(B91) is arbitrary, e.g., — q 2 = —q\ or — q 2 = —q\ [the results 
( |B9c| ) (|B9l| ) are independent of (— q 2 ) ]. In the simpler cases when only one fixed scale appears in the integrals, say 
qi = and q 2 = q, we obtain from the above results immediately 



1(0, q) 

i(o,qY 

I{0,qY ltl2 
1(0, q )W* 



1 f 1 dx 3 



(47r) 2 q 2 J v x 3 



1 



2(4tt) 2 



( V - +0) 



A 2 



-q 



1 



3(4tt) 



(g^q^ + g^ q m + jW"'^) 



hi 



A 2 



gMi qf*3 



(Blla) 
(Bllb) 
(Bile) 

(Blld) 
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As an application, using the results (B9), (Bll), we obtain the integrals on the right-hand side of DSE 



(A) 



(A) 



d A ki 1 
(2tt) 4 fc 2 (fc - kx 



r^^ 1 5/iiM2^i 2 — I( k)^ 2 kij 



2(4tt) : 



k z 



hi 



A 2 

-k 2 



d 4 fci 



k^ 



(k ki)^ (k fci)^ 



(2tt) 4 fcf(fc-fei) 2 (fc-fcO 2 



(A) Ai fc^gjjfc + gi)' 



(4tt) 2 

i 

4(4tt) : 
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A 2 
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(2tt) 4 
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(B12) 



(B13) 
(B14) 



The step from Eq. ( |bT| ) to Eq. dB14j ) follows by using the identity ( ]B2| ). 

Combining (jBl|) and flBll} ), we obtain the result for the regularized expression on the right-hand side of the DSE 



i 4Tra s (K)N c — 



( A > d 4 fc x 



(2tt) 4 fc 2 (fc-fci 



rfc" 1 



.w 2 + (£ (A) -i) 
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(fc-fci) 2 



-a s (JT) iV c 

5V ; 16tt 



(3-C (A) )ln 



A 2 
-fc 2 



(B15) 



APPENDIX C: EVALUATION OF THE VARIATIONAL DERIVATIVE OF c 



(3) 



In this Appendix, we evaluate, i.e., reduce to the form of the type (|56|), the last term in the quadratic- in-sources 
antighost DSE ( 50t ) . This term can be first rewritten as an integral over two momenta, by using the expression ( 46c ) 



/jf„c 3 (3)\ /-(A) j4i j4l 



I 6J& J 
5A b Uk 2 )W 



G c ((h+k) 2 ) 



(2tt) 4 (h + fc) 2 G A {kl)G c {{kx +k- k 2 ) 2 ) 



(h + k)^ 



SJ&(ki 



(CI) 



We note that (Sc^ /5J)(k) here means the Fourier transform of the x- function (Sc^ {x)/5J(x)). We then insert for 
the three terms in integral ( |CT| ) the ex pres sions f or c ^ 1 ' j _£_^ , A^- 1 ', in terms of the dressing functions Gx and 
the sources, as given in relations (46a)-( pffb| ) and (47a)-(47b). This results in the corresponding three integrals over 
two momenta, quadratic in the sources (~ Jp), and containing the dressing functions Ga and G c and the gluon 
propagators P^ lfij 
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y 1 >J 3 1 fc 2 J (2tt) 8 fc 2 (fci + fc) 2 fc 2 2 (fc 2 -fc) 2 GA((fc 2 -fc) 2 ) 

xfc^(fc! + k)^P^ 1 {k l )k^P^{k2)J^{k 2 )p d3 {k - k 2 ) 

_ i (K]f ac 2 c 3f c 3b2b3fb3 c 2d 3 f iA) d A hd 4 k 2 G A (fc 2 )G c ((fc 1+ fc) 2 )G^(fc 2 ) 

9 [ >J J 1 J (2tt) s k 2 {k 1 + k) 2 k 2 {k 2 -k) 2 {k 1 + k-k 2 ) 2 GA{{k2-k) 2 ) 

xfc^(fci + ky-p^Mih + k- k^p^MJ^ik^ik - fc 2 ) 

^ d 4 hd 4 k 2 G A (fc 2 )G c ((fc 1 + fc) 2 )G A ((fc 2 -fc 1 ) 2 ) 



xk^ih + k^p^M 



(27r) 8 fc 2 (fc! + fc) 2 fc 2 (fc 2 - fc!) 2 (fc! + fc - fc 2 ) 2 G A ((fci + fc - fc 2 ) 2 ) 
(2fei - fc 2 )^P A14Ml (fci) J P' I4 15 (fc 2 - fci) + (2fc 2 - kiY'P^ik^P^fa - fci) 
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The cubic terms in the structure constants can be reduced due to the following identities: 

/ aC2C3 / C3C2fc3 = -NcSab 



1 



jac 2 c 3 jc 3 b 2 b 3 jb 3 c 2 d 3 _ _ ^_ J\f ja 



b 3 , 
b2<l 3 



jac 2 c 3 jc 3 b 2 b 3 jb 2 c 2 c 3 jac 3 b 3 



(C3) 
(C4) 



Ident ity (|C3| ) is well known, and identities (C4) are direct consequences of the Jacobi identity. We use identities (C3) 
and ( |C4j ) in expressions (C2); further, we rename in the second double integral b 2 l— ► d 2 and [i 2 «-> /13; in the third 
double integral, C3 1— > <i 2 , ^3 l— ^ ^3, M4 A*5 an d /C2(new) = &2 — fci- This allows us to write expressions flC2| ) in a 
somewhat more compact form 

/J c c 3 (3)\ 

( W J (*) = 



( A ) d 4 fc lC ? 4 fc 2 G A (kl)G c ((h+k) 2 )G A (k%) f 1 fc^P M1M2 (fc 1 )(A; 1 + fc)^fc^P M 3 M4 (fc 2 ) 
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+ 2fc 2 (fc 1 +A : )^i(fc 2+ fc 1 ) 2 (fc2-fc) 2 r 2j ^ (fci) ^ ^ (fe2) 
+(2k 2 + hr-p^ik^p^ih) - (2k, + wp^MP^fo) 



fc 2 (fci + k) 2 k 2 (k 2 - kf 



x (f ad2d3 J^(k 2 )p^(k-k 2 )) . (C5) 



The next step is to perform in Eq. ( |C5|) the integration over fci. Since we work in the one-loop approximation, we need 
to obtain only the leading order term of expression (C5), i.e., the terms ~ g 3 ~ o?J 2 (note that J ~ g^ 1 ), because the 
one-loop expressions for the dressing func tions predict only the term ~ o?J 2 in Eq. ( |56| ) for the first two terms of the 
quadratic-in-sources antighost DSE ( |50b] ). Therefore, we replace in Eq. ( |Cq ) all the dressing functions Gx(k 2 ) 1— > 1. 
Nonetheless, the integration over fci is still involved, because of the rather complicated momentum dependence of the 
gluon propagators P MiMj (k s ), cf. Eq. (ff8|). However, in the Feynman gauge, £^ A ) = 1, and therefore = 1 + 0(a s ) 
by relation (|23|). So we can set £W 1— > 1 in the Feynman gauge propagators in Eq. (|C5|) at the considered order, i.e., 



Pp-itijiks) i-> Thus, we obtain 
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x (/ ad2d3 J^(fc2)p d3 (fc - fc 2 )) + 0(a^ /2 ) ■ 



(C6) 



We should first perform the integration over fci, which is complicated. Here we only take into account the cutoff- 
depen d ent t erms [~ln(A 2 /(— q 2 ))]. Then the integration of each term over fci becomes less difficult, using the results 
( B9ap (B9c) [the terms ~ln(A 2 /(— q 2 )) there ]. The integrals needed for the first, second, and third term of Eq. (Cfi) 



17 



are respectively 
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Using the results of k\ -integration (C7), we can rewrite Eq. ( |C6| ) as 
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(27r) 4 fc2(fc-A: 2 ) 2 M (4tt) 2 
(f^Jl(k 2 )p^(k-k 2 )) +0(a 5 J 2 ) 



(C7a) 



(C7b) 



(C7c) 



(C8) 



The sum in the first brackets is one. This is the result for the third term of the quadratic-in-sources antighost DSE 
(50t), in the Feynman gauge and at the leading (~ ct'J 2 ) level. 
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